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The real time evolution of field condensates is solved for 
small and large field amplitudes in scalar theories. For small 
amplitudes, the quantum equations of motion for the conden- 
sate can be linearized and solved by Laplace transform. The 
late time evolution turns to be determined by the singularities 
in the complex plane (one-particle poles, two- and multi- par- 
ticle cuts, Landau cuts for non-zero initial temperature). In 
hot scalar electrodynamics, we solve the real time evolution 
of field condensates with soft length scales ~ k~ l > (eT)~ . 
We rederive the HTL effective action using the techniques of 
non-equilibrium field theory for small amplitude condensates. 
We find that transverse gauge invariant condensates relax as 
l/t 2 and longitudinal condensates associated with plasmon 
(charge density) excitations relax with l/[t log 2 t] behavior 
to asymptotic amplitudes that are determined by the quasi- 
particle poles. 

In the nonlinear regime (for large initial energy densi- 
ties) we analyze the dynamics of dissipation and relaxation 
in the unbroken symmetry phase of scalar theory after linear 
unstabilities (parametric) are shut-off by the quantum back- 
reaction. A new time scale emerges that separates the linear 
from the non-linear regimes. This scale is non-perturbative 
in the coupling and initial amplitude. The non-perturbative 
evolution is studied within the context of the O(N) vector 
model in the large N limit. A combination of numerical anal- 
ysis and the implementation of a dynamical renormalization 
group resummation via multi-time scale analysis reveals the 
presence of unstable bands in the nonlinear regime. These are 
associated with power law growth of quantum fluctuations, 
that result in power law relaxation and dissipation with non- 
universal and non-perturbative dynamical anomalous 
exponents. 

I. INTRODUCTION 

The dynamics of relaxation of inhomogeneous field 
configurations out of equilibrium is an important prob- 
lem and a common theme in cosmology and high energy 
physics. 

In high energy physics the experimental possibility of 
studying the chiral and quark-gluon phase transition with 
high luminosity hadron colliders and upcoming heavy-ion 



colliders makes imperative the understanding of relax- 
ation and transport processes in extreme environments. 

In early cosmology, the energy and time scales during 
the inflationary stage of the universe calls for an out of 
equilibrium quantum field treatment. Moreover, the high 
energy densities involved (~ 1/g ~ 10 12 ) make necessary 
the use of non-perturbative approaches as the large N 
and Hartree methods 

There is currently a great deal of interest in under- 
standing non-perturbative real time dynamics in gauge 
theories at high temperature, both within the realm of 
heavy ion collisions and the study of the quark gluon 
plasma M- ||, as well as the possibility for anomalous 
baryon number violation in the electroweak theory. In 
both situations the dynamics of soft gauge fields with 
typical length scales > (gT)^ 1 is non-perturbative. 

Their treatment requires a resummation scheme where 
one can consistently integrate out the hard scales associ- 
ated with momenta w T to obtain an effective theory for 
the soft scales. This is the program of resummation of 
hard thermal loops [||. Physically, the hard scale repre- 
sents the typical energy of a particle in the plasma while 
the soft scale is associated with collective excitations [|| . 

Besides gauge theories, scalar models as the 0(4) sigma 
model are relevant to study the non-equilibrium relax- 
ation of an inhomogeneous initial configuration due to 
quantum and thermal fluctuations. Collective states of 
pions as those in connection with DCC (Disordered Chi- 
ral Condensates) can be treated in such a way p3| , p^[ |. 

The main task to solve in this context is the Cauchy 
problem in quantum field theory. That is, given the ini- 
tial quantum description at time t = 0, to find it for 
arbitrary time t > 0. By quantum description we mean 
the initial Fock state in case of pure states and the den- 
sity matrix when we have a mixed quantum description. 

The dynamics is determined as usual by the Fock space 
evolution equations 

91* > 

i — - — = H hJj > pure states 
dt 1 

i— = [p, H] mixed states . (1) 
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Here H is the second quantized field theoretic Hamilto- 
nian of the system and p is the density matrix in the Fock 
space. 

The expectation value of any operator O is given by 



(O) = Tr\p(t)C] . 



(2) 



For thermal initial conditions with an initial temperature 
given by 1//3, we choose the density matrix at t = to 
be 



Pi = e-* H < 



(3) 



where Hi is a free (quadratic) Hamiltonian. Notice that 
Hi is not the Hamiltonian of the system for t > 0. The 
system thus evolves out of equilibrium. The above ex- 
pectation value can be rewritten easily as a functional 
integral defined on a complex-time contour. The contour 
has two branches running forward and backward in time 
and a third leg along the imaginary axis stretching to 
t = —i/3. This is the standard Schwinger-Keldysh closed 
time path formulation of non-equilibrium field theory (see 
ref. |^,||,|l6| for details). Fields defined on the forward and 
backward time contours are accompanied with (+) and 
(— ) superscripts respectively and are to be treated inde- 
pendently. The expectation value of any string of field 
operators may be obtained by introducing independent 
sources on the forward and backward time contours and 
taking functional derivatives of the generating functional 
with respect to these sources. The imaginary time leg 
of the complex time contour does not contribute to the 
dynamics. 



II. THE AMPLITUDE EXPANSION IN SCALAR 
AND GAUGE THEORY 

Let us consider a quantum scalar field <j>(x,t) coupled 
to itself and eventually to other fields. The order param- 
eter is then defined by 



(<l>{x,t))=Tr\p{t)<Kx,t)] 



(4) 



For small values of the order parameter, we can linearize 
the evolution equations for the field <fi(x, t) with the result 

4>{x, t) - V' 2 <p{x, i) + M 2 (T) cf>(x, t) 



+ I d 3 x' I dt! T, ret (x- x',t-t') <j){x',t') (5) 



Here 'S ret (x — x', t—t') stands for the retarded self-energy 
of the field <fi(x, t) and for the renormalized mass includ- 
ing temperature corrections. Notice that eq.(||) respects 
causality since the field at time t only depends on the 
field values at earlier times t' < t. 

To be more precise, small values of the order parame- 
ter for a self-coupled scalar field of mass M means that 
y/\<j) « M where A is the coupling constant. 



In eq.(H), the self-energy kernel T, ret (x — x',t — t') can 
be computed in perturbation theory to any finite order. 
Notice that the self-energy contains information about 
the initial state through the field propagators. Namely, 
perturbative calculations are carried out with the non- 
equilibrium scalar propagators 

f d 3 k 

{$W(x,t)& b \x,t')) = -i I Gf (*,*') e -*-(*-« ) 



(2tt) 3 

where (a, b) e {+, — }. 

G+ + (t,t') = G>(t,t')e(t-t , ) + G<(t,t')Q(t'-t) , 
G^-{t,t') = G>{t,t')9(t' -t) + G<(t,t r )e(t-t') , 

G^(t,t')=-G< > \t,t'), 

G>(t,t r ) = -1- f(l + n fc ) e~ iu fc(t - t ' ) + n k e *"*<*-0' 

G k (*.*') = TT- W e-^-^ + (l + n fe ) e^C" 4 ') 
Zuj k L 



= \Jk 2 +■ 



,0u k _ I 



Here (3=1 /Ti and Tj is the initial temperature as defined 
above [eq.(|)]. 

The explicit expression and the properties of T, ret (x — 
x',t — t') obviously depend on the fields coupled to 
<fi(x,t). Several cases with scalar and fermionic fields 
were discussed in ref. Jl(],[l6| . 

Let us consider now scalar quantum electrodynamics 
(SQED). In fact, our ultimate goal is to understand re- 
laxational processes associated with off-shell effects, such 
as Landau damping, in a Non-Abelian gauge theory [El. 
What we do is to treat the same problem in the con- 
text of scalar quantum electrodynamics (SQED) model 
[Q. To leading order, we expect that this should be a 
good analogue of what happens in the Non-Abelian case; 
it also has the advantage that it is simpler to deal with, 
and as we will see below, it can be cast from the outset 
in terms of gauge invariant variables. This will elimi- 
nate any ambiguities associated with the usual problem 
of gauge dependence of off-shell quantities. 

We will start with some inhomogeneous field configu- 
ration which is excited in the SQED plasma at t = 0. 
What we want to do is to follow the time development 
of this configuration as it interacts with the hard modes 
in the plasma, and in particular, we want to know to 
determine how it relaxes in time. 

We use a gauge invariant formulation of SQED. The 
SQED Lagrangian is given by 



C = -m 2 |$| 2 



(6) 



We begin with the identification of the constraints as- 
sociated with gauge invariance. The Abelian gauge the- 
ory has two first class constraints, namely Gauss's law 
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and vanishing canonical momentum for Ao- What we 
will do is project the theory directly onto the physical 
Hilbert space, defined as usual as the set of states anni- 
hilated by the constraints. First, we obtain gauge invari- 
ant observables that commute with the first class con- 
straints and write the Hamiltonian in terms of these. All 
of the matrix elements between gauge invariant states 
(annihilated by first class constraints) are the same as 
those that would be obtained by fixing Coulomb gauge 
V ■ A = 0. The Hamiltonian in the physical subspace can 
be written solely in terms of transverse components and 
includes the instantaneous Coulomb interaction. The 
instantaneous Coulomb interaction can be traded for a 
gauge invariant Lagrange multiplier field Aq(x, t) (a non- 
propagating field whose canonical momentum is absent 
from the Hamiltonian) linearly coupled to the charge den- 
sity p(x, t) and obeying the algebraic equation of motion 
V 2 Ao(x,t) = p(x,t). This leads to the following La- 
grangian density: 

C = 9 Al $ t 9^$ + 2 d » A T ■ d ^ A T - eA T ■ Jt 
- e 2 A T ■ A T + i (VA ) 2 

+e 2 Al $ f $ - ieA a ( $$ f - , 



3t 



(7) 



where At is the transverse component of the gauge field. 

In order to provide an initial value problem for study- 
ing the relaxational dynamics of charge density fluctu- 
ations we introduce an external source 3l(x,€) linearly 
coupled to Aq and study the linear response to this per- 
turbation. Furthermore, it is convenient to introduce ex- 
ternal sources coupled to the transverse gauge fields to 
study the linear response of transverse gauge field con- 
figurations. These external fields could in principle play 
the role of a semiclassical configuration coupled to small 
perturbations in a linearized approximation. We thus 
take, 

C^C- J L {x,t)A (x,t)~ J T (x,t) ■ A T (x,t) . (8) 

The relaxational dynamics of our initial inhomoge- 
neous configurations is clearly an out of equilibrium pro- 
cess, and needs to be treated by an appropriate formalism 

The non-equilibrium SQED Lagrangian is given by 

£ noneq = £[A+,<S>+,&+,A+] -£[At,*-,&-,Ao] . 

Perturbative calculations are carried out with the scalar 
Green's functions given above in eq.(^) and the non- 
equilibrium photon propagators: 



(A^(x,t)A^(xX)) = ~ij 7^3 0#(M,O 



|(&) 



d 3 k 



ik-(x—x / ) 



g± + (k;t,t') = Vimm{t,t')e(t - o + g£(t,t')Q(t' - *)] 

Gif(k;t,t') = VijikMfrt'Mt' -t) + G£{tj)e(t - t')] 
gt j *(k;t,t') = -V ij (k)Gk i>) (t>t') 



N k = 



I 

2k 

i 

2k 



(1 + N^e-^-t"* + JV k e <fc <*- t ') 
N k e-' 1 ^-*') + (1 + N k ) e jfc (*-*') 
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e l3k _ 1 ■ 

Here Vij(k) is the transverse projection operator: 



Vijik) 



(9) 



With these tools we are ready to begin our analysis of 
non-equilibrium SQED. 

The field evolution equations in SQED for the trans- 
verse gauge field take then the form 

Aj(x,t)-V 2 Af(x,t) 

+ J dV J dt' 'E ret (x-x' ) t-t') i j Aj(x',t') = (10) 

where 1 < i, j < 3 and Yj ret (x — x',t — t')ij stands for 
the transverse photon vacuumm polarization. After cal- 
culation we find , 

d 3 p 



r ^T 2 
(1 + n p + n p+k ) 



2e* 

(2ir) 3 uj p uj k+p 



S 2 + (uJk+p + ^p) 2 



(11) 



+(n p - nk+p) 



P S 2 + (Uk+p - iOpf 



PTi PTj 



In the hard thermal loop limit (HTL) T » k, T » m 
this yields, 



e 2 T 2 f 

S fc (s)y = — — \ dQ J)Ti J>Tj 



s + ip ■ k 



(12) 



We see that the magnetic mass identically vanishes, as 
expected. Introducing the lightlike vector £ = (1, — p) 
and recognizing that the s-variable may be replaced by 
d t in the time domain, we obtained the HTL effective 
Lagrangian for the transverse electromagnetic field [Q, 



SCt 



e 2 T 2 fdfl A , 



d 



e-a 



p Tj A Tj (x,t) 



Therefore, the hard thermal loops are resummed by this 
method in a very straightforward way. 

We got analogously for the longitudinal part of the 
gauge field @, 

— A 2 (x,t)+ J —A (x,t)-J^A Q (x,t) 

which is the known result and displays explicitly the De- 
bye screening mass of charge density fluctuations. 
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III. EXACT FIELD TIME EVOLUTION BY 
LAPLACE TRANSFORM 

The field evolution equations (H) and ([h]) in the am- 
plitude approximation can be solved exactly by Fourier 
transform in space and Laplace transform in time 
LOUltJl . For the scalar case we define, 



d(k,t) 



d 3 3 



-ikx 



;2tt)3 

50 

dt e- st 



<p(x,t) , <j> k (s) 
d 3 x 



f 

Jo 



dte- st 5{k,t) 



(2tt)3 



e ip x E(jc,t) 



Notice that the Laplace transform variable s corresponds 
to i times the frequency u>. The evolution equations (|5|) 
can now be solved exactly, 



4>k{s) 



s S(k,0) + S(k,0) 
+ P + M 2 (T) + E fc (s) ' 



(14) 



We have now to perform the inverse Fourier and 
Laplace transforms: 



ix,t) 



4>k(s) 



ds 
2iri 



(15) 



where T is a contour running from —ioo to +zoo par- 
alell to the imaginary s axis and with a small but nonzero 
real part. As it is obvious the singularities of cq.(|l4|) are 
crucial for the behaviour of 4>{x,t). The zeroes of the 
denominator are specially important. 

In order to evaluate ( |l5| ) we deform the contour T 
to the left. Pole contributions are explicitly computed 
by the residue theorem. Cut contributions become real 
integrals of the corresponding cut discontinuities (see 

§!§)• 

In SQED we analogously obtain for the tranverse gauge 
field, 



ds sAf k (0)+Aj k (0)+J Ti (k,s) 



2ni 



fc 2 + ££( s ) 



(16) 



Here Ej[ (s) stands for the transverse vacuumm polariza- 
tion. 

There are two kinds of cuts in £j£(s) : 

• the two particle production cuts from s = ±i(m + 
Wk) to ±ioo where m is the mass of the scalar. 

• the thermal (Landau) cut running from — ik to +ik 

It must be stressed that the Landau cut discontinuity is 
of order T 2 whereas the production cut discontinuity is of 
order InT. Hence, only the Landau cut will be consider 
in the HTL limit. 



The explicit form for Y, k (iuj + + ) follows from eq.(ll) 
In the HTL limit we find (see for example 0) 

e 2 T 2 







u 2 ^ 


1 In 


k + u> 












k 2 






k — uj 





.e 2 T 2 ircu , 
- % —k ' ' 



k 2 



e(fc 2 - w 2 ) = £*H* + *4Mi 



In addition to the branch cut singularities, the retarded 
propagator has isolated poles corresponding to the quasi- 
particle excitations which can propagate in the plasma. 
In the HTL limit they are at s = ±iojp ~ ±*ff-- 

After deforming the contour T we obtain |7j 

A T ,(k, t) = A^I'lk, t) + A'f>(k, t) , (17) 

contributions from the quasiparticle poles add up 
to give a purely oscillatory behavior in time. The residues 
at the poles give rise to a T-dependent wave function 
renormalization. We thus obtain, 



A p T f{k,t) = Z t [T] 



A Tl (k , 0) cos u P t + A Tl (fe, 0) 

n -1 



sinujpt 

Up 



du 2 



(18) 



cj—Ljp^eT/3 



Here Z'(T) is the temperature dependent wave func- 
tion renormalization defined on-shell at the quasiparticle 
pole whose leading HTL contribution is obtained from 
the self-energy (0). The cut contribution is given by 



(19) 



A%?(k,t) = - [ dwY 1 t I (iu: + + ) 
* Jo 

A T i(k,0)uj cos(wt) +A T i{k,0) sin(wt) 



Evaluating Eqs.@,@ and 
important sum rule, 



at t = we obtain an 



2 f k 
Z t [T} + - / duj 

I" JO 



[ W 2-fc8-E*H«] 



1 (20) 



The integral over the cut (|l^) cannot be evaluated in 
closed form but its long time asymptotics is dominated 
by the end-point contribution form uj — k. We find in 
the long time limit t >> 1/k 



\ cut 



(k,t) 



12 



+ A Tl (k,0) 



3m(kt) 1 
~kl?~ J 



A Ti (k,0) 

oil 



cos(fci) 



(21) 



We want to stress that the real time dynamics of 
the condensate is completely determined by the ana- 
lytic structure of the retarded propagator and the global 
structure of the spectral density in the s-plane. 

The second important point to note is that the long 
time behaviour is a power law ~ t~ 2 (times oscillations) 
and not an exponential decay. 

This means that Landau damping effects cannot be 
reproduced by phenomenological 'viscous' terms of the 
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type ~ r 4z neither at long nor at short times. The fail- 
ure of such a phenomcnologically motivated ansatz was 
already noticed at zero temperature in different contexts 
in ref. j| . We stress that such a description not only fails 
to reproduce the power law behaviour but in fact ignores 
all the non-local physics of Landau damping which is so 
clearly encoded in the Hard Thermal Loop kernels. 

A very similar analysis yields the longitudinal compo- 
nent of the gauge field in the case of an impulsive source 
Ji,{x,t) — 5 3 (x)S(t). The case of a more complicated 
source can be obtained by convolution. In this case we 
obtain 



Ao(k, s) 



where in the HTL limit 



1 



k 2 + T, l (s) 



£'(*) 



S 2 T 2 



6 ik 



In 



is + k 



(22) 



(23) 



Thus the longitudinal self-energy along the imaginary 
axis in the s-plane, when approaching from the right is 
obtained as before to be: 



S ; (^ + + ) 



LU 

2fc ln 



k + U! 



G 



-y&(k 2 -uj 2 



(24) 



The location of the longitudinal quasiparticle or plasmon 
poles is given again for soft external momenta (k << eT) 
by s = ±iujo ~ iieT/3. The real time dependence of the 
longitudinal condensate is then found by inverting the 
transform through 



A>(M) 



e-fzoo 



ds 
2ttI 



A (k,s) 



(25) 



where the contour is to the right of all the singularities 
as in Eq. ([l6|) . As in the transverse case, at high tem- 
peratures the singularity structure is dominated by the 
discontinuity across the cut that runs from — ik to +ik 
corresponding to Landau damping, and the two plasmon 
poles at ±iluq. Deforming the contour to pick up the cut 
and pole contributions we get 



A (k,t) 



\pole 



(M) 



A CUt 



(M), 



(26) 



where 



A p ole (k,t) = -Z l [T] 



sm(ujQt) 



Z l [T] = 



dE l (iu) 



duj 2 



for the longitudinal component there is no equivalent 
sum rule because the field Aq(x) is a non- propagating 
Lagrange multiplier. 

The long time, asymptotic behaviour of the longitu- 
dinal condensate is dominated by the end-points of the 
integral (^§[). The end-point ui — yields contributions 
that vanish for long t faster than any negative power of 
t, as it was the case for the transverse part ([H)]) . We find 
that the long time asymptotics for t >> 1/fc is dominated 
by the end-point contribution at lu = k 



Afr t (k,t) = a c a f v (k,t) 



12 



5 2 T 2 



1 + - 



c dx 



cosffci + a(x, t)] 



log 4 



k/. 



5tt 2 
2 



a(x, t) = arctan 



^+log 2 f 



exp[2 



log 2 I 
6fc 2 , 



9tt± 
16 



e 2 T 2J 



The integral in Eq. ([280 cann °t be expressed in terms of 
elementary functions and it is related to the v(x) function 
ppf . For large t, one can derive an asymptotic expansion 
in inverse powers of log[/rf/c] by integrating Eq.(p8|) by 
parts: 



P(kt) = — arctan 



12 cos(fci + /3(kt)) 



e 2 T 2 t W 2 M + 2f 

° c 4 

7T log kt 

log 2 kt-4 ' 



l + O 



i 



log* 



(29) 



This expansion is not very good quantitatively unless t 
is very large. For example, for kt = 500, a™y(k,t) given 
by Eq.(|2lj|) approximates A^jk, t) up to 0.1%, whereas, 
the dominant term in Eq. ( p9[ ) is about 30% smaller than 

<fy(k,t) f. 

To summarize, we gather the final results for the 
asymptotic real-time evolution of the transverse and lon- 
gitudinal condensates in the linear approximation (a) 
Transverse (no external source): 



A Ti (k,t) = A Ti (k,0) 
+A Ti (k,0) 



Z l [T] cos Lu p t- 



12 cos kt 



e 2 T 2 t 2 



shuV 



12 sin kt 



e 2 T 2 kt 2 



O 



and 



AT(M) = — 

7T 



diu- 



E^w)/ sm(ut) 



The sum rule ( p0[ ) implies that the coherent field con- 
figuration relaxes to an asymptotic amplitude which is 
smaller than the initial, and the ratio of the final to the 
initial amplitude is completely determined by the ther- 
mal wave function renormalization. 

:T / 3 (b) Longitudinal (impulsive external source): 

(27) 

, sm(w i) 



A (k,t) = -Z l [T] 



[* 2 + ^)fl] 



"asy 



(k,t) 



l + O 



Unlike the transverse components for which the sum rule 
is a consequence of the canonical commutation relations, 



The transverse and longitudinal wave function renormal- 
izations Z*[T\ , Z l [T] are given by Eqs. aaid@ re- 
spectively and a™t(k,t) is given by Eqs. (|2Sj)-(|29|) while 



(28) 
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the positions of the poles oj p , u>q can be obtained by the 
finding the zeroes of the denominators using eqs.(|17|) and 
respectively. 

In summary, we find that the long time dynamics is 
dominated by the Landau damping thresholds at lo = 
±fc, not by the u> « region of the spectral density. The 
early time dynamics is determined by moments of the 
total spectral density. 

While an understanding of the real-time relaxation of 
non-equilibrium, inhomogeneous field configurations is of 
fundamental importance in the physics of relaxation in 
the QGP ||, such a calculation also has phenomenolog- 
ical implications for sphaleron induced B-violating pro- 
cesses. 

One of our main focus is to assess in detail the real 
time non-equilibrium dynamics of soft excitations in the 
plasma with particular attention to a critical analysis 
of the long-standing belief that the small frequency re- 
gion of the spectral function dominates the long time 
relaxational dynamics. We find, to the contrary that the 
Landau damping thresholds at ui = ±fc determine the 
long- time dynamics and that the early time dynamics 
is sensitive to several moments of the total spectral den- 
sity. In ref. ]?J we concentrated on the case of scalar 
electrodynamics (SQED) since this theory has the same 
HTL structure (to lowest order) as the non-abelian case 
jn],[l2| . Most of our results can therefore, be taken over 
to the non-Abelian case with little or no changes at least 
to the leading order in HTL. 

When finite temperature infrared divergences are 
present, like for scalar condensates coupled with gauge 
fields, the amplitude approximation combined with the 
resummation of the infrared divergences by dynamical 
renormalization group leads to anomalous logarithmic re- 
laxation and logarithmic phases in time ]2l| . 

IV. NON-LINEAR QUANTUM EVOLUTION IN 
SCALAR FIELD THEORIES 

Our goal is to deal with the out of equilibrium evolu- 
tion for large energy densities in field theory. That is, a 
large number of particles per volume m -3 , where m is the 
typical mass scale in the theory. The most familiar tech- 
niques of field theory, based on the S-matrix formulation 
of transition amplitudes and perturbation theory apply 
in the opposite limit of low energy density and since they 
only provide information on in — > out matrix elements, 
are unsuitable for calculations of time dependent expec- 
tation values. 

Similar tools are also necessary to describe consistently 
the dynamical processes in the Early Universe [QJ|]. In 
particular it has been recognized that novel phenom- 
ena associated with parametric amplification of quan- 
tum fluctuations can play an important role in the pro- 
cess of reheating and thermalization |l],|3|,[l7]] . It must 
be noticed that the dynamics in cosmological spacctimes 



is dramatically different to the dynamics in Minkowski 
spacetime. Both in fixed FRW and de Sitter back- 
grounds and in a dynamical geometry the dynamical 
evolution is qualitatively and quantitatively different to 
the Minkowski case considered in the present paper. 

Our program to study the dynamical aspects of re- 
laxation out of equilibrium both in the linear and non- 
linear regime has revealed new features of relaxation in 
the collisionless regime in scalar field theories ■ Re- 
cent investigation of scalar field theories in the non- linear 
regime, including self-consistently the effects of quantum 
backreaction in an energy conserving and renormalizablc 
framework have pointed out to a wealth of interesting 
non-perturbative phenomena both in the broken and un- 
broken symmetry phases ||. These new phenomena are 
a consequence of the non-equilibrium evolution of an ini- 
tial state of large energy density which results in copi- 
ous particle production leading to non-thermal and non- 
perturbative distribution of particles. Our studies have 
focused on the situation in which the amplitude of the ex- 
pectation value of the scalar field is non-perturbatively 
large, A w \/A < $ > /m ~ 0{\) (m is the mass of 
the scalar field and A the self-coupling) and most of the 
energy of the initial state is stored in the 'zero mode', 
i.e. the (translational invariant) expectation value of the 
scalar field $. Under these circumstances the initial en- 
ergy density £ « m 4 /A. During the dynamical evolution 
the energy initially stored in one (or few) modes of the 
field is transferred to other modes resulting in copious 
particle production initially cither by parametric ampli- 
fication of quantum fluctuations in the unbroken symme- 
try phase, or spinodal instabilities in the broken symme- 
try phase. Furthermore, we have considered in ref. |l4| ] 
the evolution of an initial state with a large number of 
quanta with a distribution around a momentum |fco| cor- 
responding to a thin spherical shell in momentum space, 
a 'tsunami'. 

In both cases the mechanism of energy dissipation and 
particle production results in a number of produced par- 
ticles per unit volume M oc m 3 / A which for weak cou- 
pling is non-perturbatively large |||l4|]. We call 'linear 
regime' to this first stage dominated either by paramet- 
ric or spinodal unstabilities. 

We recognized || a new dynamical time scale t\ where 
the linear regime ends. By the time t\ the effects of the 
quantum fluctuation on the dynamical evolution become 
of the same order as the classical contribution given by 
the evolution of the expectation value of the field. The 
'non-linear' regime starts by the time t\. In the case 
of broken symmetry, this time scale corresponds to the 
spinodal scale at which the backreaction of quantum fluc- 
tuations shut-off the spinodal instabilities. At this scale 
non-perturbative physics sets in and the non-linearities 
of the full quantum theory determine the evolution. This 
time scale t±, which we call the non-linear time, is a non- 
universal feature of the dynamics and depends strongly 
on the initial state and non-perturbatively on the cou- 
pling, as ti cx log[A _1 ] for weak coupling ||. 



G 



There are two very important parameters that influ- 
ence the quantum dynamics: the strength of the coupling 
constant A and the initial energy density in units of the 
scalar field mass m. If in the initial state most of the 
energy is stored in the zero mode, the energy density is 
determined by the amplitude A of the order parameter. 
The value of this field amplitude determines the regime 
of applicability of perturbation theory methods. That is, 
the amplitude approximation considered in sees. II and 
III fails to catch important aspects of the physics when 
A is not small. 

We consider now the non-perturbative regime in 
which A = VA < $ > /m ~ 0(1). It is important 
to point out that for large field amplitude, even for very 
weakly coupled theories non-linear effects will be impor- 
tant and must be treated non-perturbatively. This is the 
case under consideration. Having recognized the non- 
perturbative nature of the problem for large amplitudes 
we must invoke a non-perturbative, consistent calcula- 
tional scheme which respects the symmetries (continuous 
global symmetries and energy- momentum conservation), 
is renormalizable and lends itself to a numerical treat- 
ment. 

We are thus led to consider the O(N) vector model 
with quartic self-interaction || and the scalar field in 
the vector representation of O(N). 

The action and Lagrangian density are given by, 



S = 



8N 



, < \ [d^(x)f - V(*(x)) 
, 2Nm 2 \ 2 N to 4 



A 



2A 



(30) 



Heisenberg equations of motion ]3|. [l5[|l9| ] given below ( p2| - 
0). 

Our choice of initial conditions on the density matrix 
is that of the vacuum for the instantaneous modes of 
the Hamiltonian at the initial time P,19|. Therefore we 
choose the initial conditions on the mode functions to 
represent positive energy particle states of the instanta- 
neous Hamiltonian at t = 0, which is taken to be the 
initial time. That is, 



A(o) = 



l 



; A(0) = -iVw k ; W k = Jk 2 + M 2 



where the mass Mq is defined below. With these bound- 
ary conditions, the mode functions /fc(0) correspond to 
positive frequency modes (particles) of the instantaneous 
quadratic Hamiltonian for oscillators of mass Mq. 

We point out that the behavior of the system depends 
mildly on the initial conditions on the mode functions as 
we have found by varying them within a wide range. In 
particular, the various types of linear and nonlinear reso- 
nances are independent of these initial conditions PJl9| . 

We introduce the following dimensionless quantities: 



= H t > 1 = T— 7 > n q = T~~7 ' ^ 2 ( T ) 



A 



2 m 



$ 2 (i) 



2 m 



[<* 2 (*)>h-<* 2 (0)>h] , (£(0) = 0) 



Here (^ 2 (t)) n stands for the renormalized composite op- 
erator [see eq.(B5J) for an explicit expression]. 



We write the field <fr as # = (cr, tt) where n represents 
the TV — 1 'pions', and choose the coupling A to remain 
fixed in the large N limit. Due to lack of space we only 
considered here the unbroken symmetry case to 2 > 0. 
The broken symmetric case (to 2 < 0) is considered in 
ref. |. 

We can decompose the field a into its expectation value 
and fluctuations x(a:,i) about it: 



(31) 



a(x,t) = $(t) VN + x(x,t) , 



with $(t) being a c-number of order one in the N — > 00 
limit and x an operator. 

We shall not rederive here the field evolution equa- 
tions for translationally invariant quantum states. In the 
leading order in the large N approximation the theory 
becomes Gaussian at the expense of a self-consistent con- 
dition ||[ij|,]l9|, this in turn entails that the Heisenberg 
field operator *ff(x,t) can be written as 



9(x,t) = 



d 3 k 1 
(2^7! 



a u f k (t) e ik -* + alf k (t) e 



where ak , a k are the canonical creation and annihilation 
operators, the mode functions fk(t) are solutions of the 



A. Evolution equations in the large N limit 

In terms of the dimensionless variables introduced 
above the renormalized equations of motion are found 
to be (see references ) : 



V + V + if + 9 V(r) S(r) = 

" d 2 



dr 2 



+ q A + l+ ry(r) 2 +.gS(r) 



<P 9 (t) = 



(32) 
(33) 



<P q (0) = 



<p q (0) = -ijn q , t?(0) = 770 , jj(0) = (34) 



Hence, M 2 {t) = 1 + ^{t) 2 + g £(t) plays the role of a 
(time dependent) effective mass squared. 

We will choose J7 g such that at t — the quantum fluc- 
tuations are in the ground state of the oscillator s at the 
initial time. Since <?£(0) = 0, we choose Mq = \/ 1 + rj 2 . 
g£(r) is given by the self-consistent condition p,|l^| 



9 E(r)=g J™ q 2 d q ^ q {r) \ 2 -i- 



2q 3 



V 2 (T)+gj:(r)] 



(35) 
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We thus see that the effective mass at time r contains 
all g-modes and the zero mode at the same time r. The 
evolution equations are then nonlinear but local in time 
in the infinite N limit. 



B. The early time evolution: parametric resonance 

In the weak coupling regime the back-reaction term 
ffS(r) is small for small g during an interval say < r < 
Ti P,|l9[. This time n will be called the nonlinear time 
and it determines the time scale when the backreaction 
effects and therefore the quantum fluctuations and non- 
linearities become important. 

For times t < Ti, 3S(t) can be neglected and then 
eq.(p2) reduces to the classical equation of motion 



fj + r] + 77' 1 = 



(36) 



The solution of this equation with the initial conditions 
can be written in terms of elliptic functions 



77(1-) = 770 cn (r\Jl + r? 2 , > 



k = 



where cn stands for the Jacobi cosine. 

Inserting this form for t](t) in cq.(|33"|) and again ne- 
glecting gS(r) yields 



+ <f + l + 77 2 cn 2 [rJl + rjlk 



dr 2 



tp q (r) = 



(37) 



This is the Lame equation for a particular value of 
the coefficients that make it solvable in terms of Jacobi 
functions IH). 

Since the coefficients of eq. ( |37| ) are periodic with period 
2oj (notice that t](t + 2lo) = —j](t)), the mode functions 
can be chosen to be quasi-periodic (Floquet type) with 
quasi-period 2oj, 



U q (T + 2Lo) = e* F ^ U q {r), 



(38) 



where the Floquet indices F(q) are independent of r. In 
the allowed zones, F(q) is real and the functions U q (r) 
are bounded. In the forbidden zones F(q) has a non-zero 
imaginary part and the amplitude of the solution either 
grows or decreases exponentially. 

It turns out in the present case that there is only one 
forbidden band for q 2 > running in the interval < 
<Z 2 <% 2 /2. 

The modes from the forbidden band < q < rj a /\/2 
dominate S(r). For < r < n, S(t) oscillates with an 
exponentially growing amplitude. This amplitude (enve- 
lope) Ti env (t) can be represented to a very good approx- 
imation by the formula [p|,[19[ 



,{r) 



N < 



(39) 



where B and N are functions of 770 given by 
B(») = 8 \/l + r$ q (1 - 4q) + 0(f) , 



NM = JL ^<i±l|h^M |1 + 0W)1 . (40) 

and the elliptic nome q can be written as a function of 
i]o as 



1 (l+Tfa 2 ) 1 / 4 - (l + r/ 2 /^/ 4 

2 (1+% 2 )V4 + (1+^2/2)1/4 



(41) 



with an error smaller than ~ 10 -7 . 

Using this estimate for the quantum fluctuations S(r), 
we estimated the nonlinear time scale t\ at which the 
back-reaction becomes comparable to the classical terms. 
Such a time is defined by gT,(n) ~ (1 + r? 2 /2), @ 



1 



n 



B(vo) 



log 



.iYM(l + % 2 /2) 




FIG. 1. The zero mode 7y(r) vs. t for the unbroken sym- 
metry case with t)q — 4, g = 10 -9 . 

The time interval from r = to r ~ T\ is when most 
of the particle production takes place. Around r ~ T\ 
the quantum fluctuation become large enough to begin 
shutting-off the growth of the modes and particle pro- 
duction slows down dramatically. This dynamical time 
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scale separates two distinct types of dynamics, for r < T\ 
the evolution of the quantum modes tp q (r) is essentially 
linear, the backreaction effects are small and particle pro- 
duction proceed via parametric amplification. For r > t\ 
the quantum backreaction effects are as important as the 
tree level term t](t) 2 and the dynamics is fully non-linear. 

The growth of the unstable modes in the forbidden 
band shows that particles are created copiously (~ 1/g 
for t ~ Ti). Initially, (r = 0) all the energy is in the clas- 
sical zero mode (expectation value). Part of this energy 
is rapidly transformed into particles through parametric 
resonance during the interval < r < t\. At the same 
time, the field expectation value decreases as is clearly 
displayed in fig. 1. 

The momentum distribution of the particles produced 
by parametric resonance follows the Floquet index and is 
peaked at q w ^t/q (1 — q) fl|l9|. This is shown in fig. 2. 




FIG. 2. Momentum distribution of the produced particles 
n(q) = q 2 Nq d (T) at the nonlinear time r = r\ for rjo = 4, 
ff = l(T 9 



non-perturbatively large amplitude ~ 1/ VA which could 
be described as a 'semiclassical condensate' in the unbro- 
ken phase [see fig. 3]. 




FIG. 3. Momentum distribution of the produced particles 
at r = 3000 for rjo = 1, g = 10 -3 . Notice the main peak at 
q — 0.067, in good agreement with the estimate Ki/t given 



m 



ref.[2]. 



The novel result that emerged combinining numeri- 
cal methods and dynamical renormalization group anal- 
ysis is the presence of non-linear resonances that lead to 
asymptotic relaxation described by non-universal power 
laws \ |] . These power laws are determined by dynamical 
anomalous exponents which depend non-perturbatively 
on the coupling. 



For late times, typically larger than 3 or 4 times ti, the 
effective mass tends asymptotically to a constant value 

2 

Ai^e = 1 + -j-- This fact implies that the modes becomes 
effectively free. 

At the same time, the field expectation value relaxes 
to zero asymptotically with a non-universal power law. 
The initial energy density which is non-perturbatively 
large goes completely into production of massive parti- 
cles. The asymptotic particle distribution is localized 
within a band determined by the initial conditions with 
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